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O , Abstract 

a" 

We study the physical consequences of two diffferent but closely re- 
lated perturbation schemes applied to the Einstein-Maxwell equations. 
In one case the starting space-time is flat while in the other case it is 
Schwarzschild. In both cases the perturbation is due to a combined elec- 
tric and magnetic dipole field. We can see, within the Einstein-Maxwell 
equations a variety of physical consequences. They range from induced 
gravitational energy-momentum loss, to a well defined spin angular mo- 
J- ' mentum with its loss and a center-of-mass with its equations of motion. 

en 

O ■ 1 Introduction 

O 

Recently, using the spin-coefficient (SC) formalismp], a perturbation scheme 
for a simple model using the Einstein-Maxwell equations was described^ GO- 
Pure electric and magnetic dipole radiation was considered as a first order per- 
turbation off a flat space background. It was shown that this in fact leads to 
a second-order perturbation in the gravitational field (the Weyl tensor) , and in 
turn a perturbed second order metric. Some interesting results were obtained 
by simply looking at the asymptotic Weyl tensor; these included the existence 
of a Bondi news function created by the dipole radiation with the accompanying 
classical Bondi gravitational and electromagnetic energy- loss. A pretty result 
was that one could identify, in the Bianchi Identities, the classical electromag- 
netic angular momentum loss [4]. 

In the present work, we apply this scheme to more complicated perturba- 
tions. We find, perturbatively, two different versions of what could be loosely 
characterized as generalized Reissner-Nordstrom space-times; that is, metrics 
with a mass and Coulomb charge, but now with electromagnetic dipole radia- 
tion. In one case (perturbations off Schwarzschild) we consider the mass to be 
zeroth order and the charge and dipole fields to be first order. In the second 



case (perturbations off flat space) we consider the mass, the charge and dipole 
field all to be first order. 

In Section II, the Maxwell equations, for both perturbation types, are first in- 
tegrated. (It does not actually matter in which background [flat or Schwarzschild[ 
this integration occurs for these two cases.) Next, (Sections III and IV) with 
the Maxwell field in the stress tensor as the source, we integrate the Bianchi 
identities obtaining the radial and non-radial behavior of the Weyl tensor. We 
then probe further the asymptotic behavior of the Weyl tensor, in particular 
looking at the angular momentum loss and Bondi energy-momentum loss the- 
orem as well as dynamical equations for the motion of the center of mass and 
charge. Different physical consequences for the two different perturbations are 
discovered. 

For completeness, in Section V the full behavior of the second order metric 
is presented. Sec. VI contains the discussion, while an appendix contains full 
expressions for the Weyl tensor and spin coefficients. 

2 The Maxwell Field 

We work in the Bondi coordinate system (u,r,(,() , where u labels the light 
cones, C, with apex on a time-like world-line, r is the affine parameter along 
the null geodesies, and £ = cot(#/2)e*^ is the complex stereographic angle la- 
beling the null geodesies on <£. Furthermore, we choose the Bondi null tetrad 
{l a , n a , m a , m a } such that the vector l a is tangent to the null geodesic congru- 
ences. This, along with the choice of l a as a gradient, fixes the spin coefficients 
K,ir,s,p, and r as [5j: K — ir = e = 0, p = p, and r = a + (3. In the 
Schwarzschild space-time, this fixes the set of SCs as[6]: 

K = 7T = £ = CT = T = !/=A = 0, (1) 
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The flat space Minkowski set of SCs are obtained simply by setting the 
Schwarzschild mass, Ms, equal to zero. In both cases it can be seen that the 
radial and non-radial Maxwell equations are the same[5l[7], allowing us to write 
down the desired solution which contains both a Coulomb charge and radiating 



electromagnetic dipole given by: 
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The three-vector D 1 is our complex dipole moment and can be written as 
the complex superposition of its real electric and magnetic 
parts: 

D*=D l E + iD* M , (3) 

while q is our Coulomb charge. At this point, D l may be regarded as an 
arbitrary function of the retarded time, u r = \/2u. In both perturbation cal- 
culations that follow, we will treat both D % and q as being first-order quantities 
in the perturbation. All calculations will be performed to second order. In 
addition we only keep the I = 0,1,2 spherical harmonics. 

Remark Other papers [U [10] have often written D l = qS, 1 , where £* is a 
complex position vector. Due to our perturbation formalism, we can make no 
such identification at this point in the calculation. 

Remark Throughout this paper, we denote differentiation with respect to 
Bondi time as: d u () — ('), while differentiation with respect to the retarded 
time, u r = \/2u 7 is given as d Ur Q — ('). Thus, we see that (') = v2('). In later 
sections, when we want to restore units where c ^ 1, we must take (') — > c _1 ('). 
The gravitational coupling constant is k = 2Gc~ 4 . 

3 Schwarzschild Background Perturbation 

We now consider the Maxwell field given by Eq. |(2|) to be a first order perturba- 
tion in the background of the Schwarzschild space-time with the spin-coefficients 
given by Eq.([T]). In what follows, we carry this information, via the stress ten- 
sor, into the Weyl tensor by integrating the SC form of the Bianchi identities. 
Looking at the asymptotic Weyl tensor components we study the physical con- 
sequences (i.e., mass, momentum, angular momentum and equations of motion 
seen at null infinity) of our model. 

3.1 The Radial and Non-radial Bianchi Identities 

We are seeking a solution of the Bianchi Identities which is driven exclusively 
by the original Schwarzschild mass (zeroth order) and electromagnetic pertur- 
bation; i.e. with no gravitational degrees of freedom. This leads to ip = 0. 



(This step has been taken by others, e.g., [21 E]). The radial Bianchi identities 
are then given by [5]: 
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The differential operator A is given by: A = d u — d r -|- r ~c _ zv zuivisu r 
Integration produces the following results. (We relegate the complicated 
terms, Ai, involving higher r-dependence to the appendix, as these are largely 
unnecessary for the calculations of particular interest in this paper) : 
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^3 = 4+^3 (7) 
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04 = ^ + M. (8) 

Here, the 05,02,03 and 04 are r - independent functions of integration, 
which can be determined from the non-radial Bianchi identities, which take the 
form [5] : 

80? = 3fc0°0°, (9) 

i>l = -80° + 2^°>°, (10) 

02 = -303 + ^202, (11) 

0° = -a^2- (12) 

The integration of Eqs.(f9])- ([12|) is relatively straight forward, basically entail- 
ing the comparison of spherical harmonic coefficients. We make considerable 
use of Clebsh-Gordon expansions [TT] of the quadratic terms in the Maxwell 
field. This process applied to Eq.((9]) yields: 

^ = 3kD i & ,, Y^ j +if, 1 k Y^. (13) 

The vector tp 1 (u r ) emerges as an unknown integration factor in much the 
same way that if) x emerged from the integration of the radial equations. Shortly, 
however, we will be able to determine 0° up to constants. From Eq. lflOJI . we 
find: 



tfS = -^-Ms + T.-H^;*' + 2kiD>"IP' tlJk )Y° k (14) 



where the I = function Y e (u r ) is a (first-order) function of integration. Con- 
sidering Eq.(|lip. we recognize that as 3 is a spin weighted s = — 1 quantity it 
has no I = harmonic contribution. This allows us to both obtain 3 as well 
as place a restriction on T e : 

0° = (V2kiD i "Di"e ijk + 2V2ki(D i 'D j "ye ijk -^"^Y^ 1 (15) 

-2y/2kqD k "'Y{? + k (\D t "D : >" - (D l &")" - hjy'D'")'] Y^}, 
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Finally, we determine both V'i and ip 4 from Eq,( [T2"|) . We recall that ip 4 is 
an s = —2 quantity and hence does not contain an I — 1 harmonic. The I = 1 
harmonic contributions to the equation thus must vanish yielding a differential 
condition on ip° (u r ); the remaining part yields the determination of ijP 4 : 
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For later use we decompose Eq. l|18|l into its real and imaginary parts using 
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V>%"' - 2V2k(D^D%y"e ijk -V2k(D%D%yei jk -2V2kqD%"' (19) 
^f" - 2y/2k[(D%l%' + n&D%)e ijk + qD%]''. (20) 

Note that the latter equation can be immediately integrated twice as 

^' = 2V2k{(D%D% + D%D%)e ijk + qD%], (21) 

where we have taken the two constants of integration to vanish while the first 
can be integrated once as: 



$ % = V2k[2(D j ^D%)" - D%'D%]e ijk - 2V2kqD»» . (22) 

Remark: Bramson [3], improperly set the ip x = 0, apparently overlooking 
this differential condition. 

Remark: There is a further equation (a reality condition) on i\)\ that plays 
a very important role. It however involves a further variable (the spin-coefficient 
er°, i.e., the Bondi shear) for its description. 



3.2 Reality Conditions 

To more easily see the physical content in our equations, we now introduce the 
c in all time derivatives, i.e., (') — > c~ 1 ('). 

We turn to the Bondi mass aspect where its reality forces further restrictions 
on our variables. To proceed, we must first know the shear i.e., the spin 
coefficient a. A more detailed discussion of the spin coefficients is given in 
Section V and the Appendix, but it can be seen [2] that the behavior of a does 
not change with the switch from flat-space to a Schwarzschild background or 
the addition of a Coulomb charge. Thus, we find[2]: 
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Using Eqs, f23|) and l|25p . the Bondi mass aspect OH] is given by 

* = ^ + 5V + \/2c- 1 ( 7 ( ( 7 )', (26) 
and is subject to the reality condition: 

* = *. (27) 

As can be seen from Eq. (|25|) . the quantity a is second order in the per- 
turbation, so that it can be neglected in the definition of the mass aspect in 
Eq. (|26| . If we expand the mass aspect in spherical harmonics, 

* = *° + &Y& + ¥ j Y^ + ■■■ , 
then from BondiOH], the I = and I — 1 terms, i.e., 

*° = -^M g + T e = -^M B (28) 

¥ = T P* (29) 

are, up to numerical factors, interpreted as the Bondi mass and linear three- 
momentum respectively. 

From Eq.([2Sl) using Eqs.JUl) and ([25]), we have that: 

* = -^Ms + T e +(^ 1 k ' + ^&''D^ ijk )Y 1 k (30) 
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By the symmetry on the I = 2 contribution, it follows that \£* J is real. From 
the I = contribution, we see that the Bondi mass is real: 

2V2G, , „ 2V2G, r - 2^26? „ r , , 

- -^2— M s + T e = — -o- M s + T £ = — -^M B . 31 

c c c 

The reality condition on the I = 1 condition is a bit more complicated. First, 
we decompose the dipole D l into its electric and magnetic parts, 

and write: f/'i = ip 1R +iipu . Then the condition for \fr* = & yields two 
relations; one on the real and one on the imaginary part of i/'i • We find for 
the vanishing of the imaginary part ip® , 

i>l? = 2V2kc- 1 qD^ + 2V2kc- 2 (D%D% + D&D%)e iJk , (32) 



which is identical to the earlier derived, Eq. ([21]) . The real part leads to an 
expression for the Bondi linear momentum 
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(33) 



(34) 



By taking the u r derivative and simplifying via Eq, l(22|) . we obtain the electro- 
magnetic momentum flux law, 
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The imaginary part of the / = 1 harmonic of ip 1 , i.e., ^p u , is often, in vacuum 
linear theory, taken as proportional to the total source angular momentum as 
viewed from infinity. This becomes modifiedpQ]] in the presence of a Maxwell 
Field as 
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Eq. ([32|) is seen as the classical law of conservation of angular momentum for 
electromagnetic dipole radiation[4\ or angular momentum flux law: 
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3.3 Null Rotations and Equations of Motion 

We now construct a transformation (a null tetrad rotation around n a ) to what 
we define as the complex center of massJlOty. 

Though this is not the place to go into a detailed explanation |13l ITOl [9j of 
the meaning of the term ' complex center of mass ', nor into the details how it 
can be calculated or found, a brief explanation is in order. 

In a given asymptotically flat space-time, the family of regular asymptot- 
ically shear-free null geodesic congruences are determined by 1. the Bondi 
asymptotic shear a (u r ,(, C) an d 2. an arbitrary choice of a complex world- 
line that 'lives' in the space of complex Poincare transformations (complex 
Minkowski space), a subgroup of the BMS group acting on 3 + . {That an asymp- 
totically shear free null geodesic congruence picks out a complex world-line in 
complex Minkowski space is the central fact in the present discussion. That 
it has led to a series of remarkable results[10l [14j is the defense [j~5] of its rel- 
evance.} A particular complex world-line can be chosen so that the asymptoti- 
cally defined center of mass and angular momentum both vanish on it. The basic 
variable used to describe the asymptotically shear-free null geodesic congruence 
is a (stereographic) angle field, L(u r , £, () on 3 + that points backwards into the 
space-time determining a past null direction at each point (u r , C, C) of 3 + . The 



angle field, for asymptotically shear-free congruences, satisfies the differential 
equation [13] 

bL + LL, u =a°. (38) 

The solution to Eq. (f3"8"]l . (accurate to our working order), is 

L(u r , C, C) = e y?i - l?v% jk Y& + • • • , (39) 

where C( u r) and v l = £ l are respectively the arbitrary complex world- line and 
its velocity. Quadratic terms and high harmonics [10] have been omitted. 

The transformation (null rotation) on 3 + oi the Bondi tetrad ( l,n,m,fh) 
to a new tetrad (l*,n*, m* , m*) (where the I* is the null tangent vector to the 
asymptotically shear-free null geodesic congruence) is given by 



I* = I + —m + -to + 0(r~ 2 ) (40) 

r r 

TO* = m + Oir^ 1 ) 
n* = n, 

Remark We have gone from the Bondi tetrad frame with null vector I that 
is twist free but has shear to a null vector I* that is asymptotically shear-free 
but now possesses twist. The twist, S, given by 

tS = -(3L+ V2LL' - SL - V2LL'). 

It vanishes, in our approximation, when the world-line £ l is real. 

This transformation induces a transformation of the asymptotic Weyl tensor 
components. In particular the Weyl component tp 1 transforms as 

V^* = ^-3LV^ + ---, (41) 

The basic physical idea is that in linear theory the I — 1 harmonic component 
of ip 1 is (usually) taken as proportional to the complex center of mass, i.e., the 
(real) center of mass + i angular momentum. Our procedure is now to chose 
the arbitrary complex world-line so that the I — 1 harmonic component of tp-j* 
vanishes. The world-line so obtained is the complex center of mass. We thus 
have to solve 

= rl>° li -3Lil%\ i + ..- (42) 

for £ s (u r ). 

Using Eqs. |(13]), JH]), and |(39]), Eq.g2]) becomes: 

= V^ + 6 -^M s e - 3 -^M s ZVe ijk . (43) 



From Eqs. (jTSJ) , 1(32]) . and (J33J) , we saw that at least up to initial conditions, the 
vector tf> 1 is a second order quantity in our perturbation framework. However 



by choosing time-independent first-order initial conditions, V'i(o)) we see from 
Eq.((43|, that 

k _ V2c 2 V^f 0) 

€ ~ 12GMs ' (44) 

i.e., £ is a constant vector and complex center of mass is at rest. 

In this case, where £ is a constant vector, the real part can be set to zero 
by a Poincare translation on 3 + so that 

2 I Ok 

k=e _^i(o L _ 
6V2GM S 

Thus, we see that when the mass of the system is considered to be zeroth 
order, the equation of motion for the center of mass is trivial: the center of 
mass simply sits on the time axis. Physically, this can be thought of in the 
following manner: if Ms is the initial Schwarzschild mass, then it is too "heavy" 
for its motion to be affected by the "small" electromagnetic perturbation given 
by Eq.|f2]) (at least to second order in the calculation). We see later that for 
small mass the situation is very different. 

3.4 Physical Interpretations 

It turns out that there are a variety of physical interpretations - some new and 
some old - to the results of this section. 

First of all when the appropriate units are inserted, Eq. (fl6|) is exactly the 
Bondi mass/energy loss equation 

M' B = -1L(D»'D% + D^D)%). (46) 

At this approximation it coincides with the classical electromagnetic dipole en- 
ergy loss. At another approximation level there would be a forth-order correc- 
tion for gravitational energy loss via the square of the Bondi news function, i.e., 
quadrupole radiation. 

The Bondi momentum loss, given by Eq. l(35|) . 

P» = ±D%D%e ijk (47) 

is just the electromagnetic momentum flux. 
More interesting is Eq.JM), (or (|2l)l ): 

$? = 2V2kqc- 1 D^ + 2V2kc- 2 (D%D% + D%,D%)ei jk . (48) 

The imaginary part of the I = 1 harmonic oiip lJ i.e., ip u , is often, in vacuum 
linear theory, taken as proportional to the total source angular momentum, J fc , 
as viewed from infinity. This becomes modified[10j in the presence of a Maxwell 
Field as 

10 



Eq. (|48j) is seen as the classical law of conservation of angular momentum] '41 
for electromagnetic dipole radiation: 

J k ' = ^( d eD 3 e + D%D^)e ijk . (50) 

In the past, the imaginary part of the world-line vector, £j, has been identified [TBI 
[TO] with the intrinsic spin associated with the asymptotic metric via the rela- 
tionship 

S* = MbcCi- (51) 

We then see, from Eq. l[44j) . that the initial value of ip x , i.e.,ipu \ is propor- 
tional to the (constant) spin 

c 3 tb 0k 

»**—&& = *■ (52) 

so that the angular momentum, J k , in Eq. l(49l) . consists of the sum of three 
terms, the intrinsic spin, S z , a complicated term involving the integral over 
quadratic derivatives of the dipole moments and an unusual term proportional 
to the derivative of the magnetic dipole moment, ^c~ 1 qD k ^ I . 

4 Minkowski Background Perturbation 

We now consider an alternative perturbation scheme, namely perturbations off 
Minkowski space-time. The gravitating mass now enters into the calculation 
as a first order quantity in the mass term of the Weyl tensor component ip 2 . 
In the previous model the mass was large; now it will be small. Most of the 
calculations of the previous section are essentially identical in this framework, 
so we can proceed quickly toward the null rotation calculations where new results 
do appear. 

4.1 Radial and Non-radial Bianchi Identities 

Once again, we set tp — so that our solution is driven purely by the per- 
turbative quantities. The radial Bianchi identities are simply those in flat 
Minkowski space-time. The results of their integration are the same as those 
obtained earlier[2j (but now with the addition of a Coulomb charge) or by setting 
M s = in Eqs.gj). 

At the leading order in r, this yields radial behavior equivalent to that found 
in the preceding section. The non-radial integrations give virtually identical 
results: 

4>i = $ "Yil + Zkc- 2 D i &"Y^ , (53) 
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^o = Te+ l^ D k» + ^o k , + ^ & ,, DJ , eijk]Y O k (54) 

V2k I ' [D l &")' D l "Di n 



Y° 



$ = \^&''D^ jk+ 2 ^(D i '&'ye ijk -c- 2 ^'yr k 1 (55) 

C V o o / 

T e = V2T' e = ^D<"S»'"^, (56) 

^° = 72fcc- 5 f(£W")"' + l(D tf &")" ~ ±(&"Di"y\ Y~f, (57) 

^ k " = V2kc- 2 i[(D i "D j ")e ijk + 2(D i, D j ")'e ijk + 2iqcD k '"}. (58) 

We have denoted the mass term of the perturbation as T e to indicate that 
it is a first-order, perturbative quantity which arises as an integrating factor in 
tp2- The absence of the Schwarzschild mass is the essential change from the 
previous section. We now obtain the reality conditions for this new calculation. 

4.2 Reality Conditions 

As ip 4 has not changed from that obtained in Section 3.1, it follows that the 
value of the spin coefficient a will remain that given in Eq. ([25j) . so the Bondi 
mass aspect is simply given as: 

* = T e + (^D k " + ^? fc/ + 2 -^D i "tf>e ijk )Y? k (59) 

+^j(D i "&")du r Y 2 %, 

and consequently, the reality conditions are little changed from those of Section 
3.2. Once again, the I = 2 condition is trivially satisfied, while the Bondi mass 
and linear momentum are given respectively by: 

Te = T e = - 2 -^-M B (60) 

* fe = - 6 4 pk = &!' + #^. - 2 4iKD%)^l (61) 

cr c Ac c' 
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The reality of the I = 1 coefficient of Eq. ([59j) again yields: 

$? = 2V2kqc- x D% - 2V2kc- 2 (D%D^ + D%D^)e ijk , (62) 

Aside from Eg. (|60|1 these are identical to those obtain in the prior section. 

4.3 Null Rotations and Equations of Motion 

It is in our attempts to understand the physical content of our equations, i.e., 
conservation laws, the definition of angular momentum and spin, and equations 
of motion for complex center of mass and charge world- lines that the present per- 
turbation scheme departs from that of the Schwarzschild background. In par- 
ticular, using the same stereographic angle field L(u, £, C) described in Eg. (|39)1 
and applying it to the null rotation of ip li given in Eq. lj42"Tl . we obtain (since 
Mb, q and £ are first order) the second order relation: 



= ^ + 6V2G M ^ fc (63) 



or decomposed as 



i>1% = ~^M B & (64) 

c 

i ok 6v2G 

&■ 

Taking the prime derivative of Eg. (|65)1 and inserting it into Eq. (|62|) we obtain 
after simplifications, again the classical conservation law of angular momentum[4], 

(aW + ^- 2 <zO' = i^ + flgDJiJei,* (66) 

rkl 



J ' — ang .mom. flux 



with the identifications 



S k = M B ctf, (67) 

J k = S k + f c D%. (68) 

This can be considered as the evolution equation for the spin, S k . 

In order to obtain the dynamical law for the (real) center of mass £ R , we 
first note, from Eq. lj64j) . that 

r 2 

MBiR - _ 672G^ lfl ' 
Then, with the use of Eq.(j22]), i.e., 
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$ %' = V2kc- 2 [2(D^D%)" - D%'DQe ljk - 2^k C - 1 qD%", (69) 

it becomes 

Mb&' = -^[(D%D^)" - \D^DZ]e l]k + J^l'", (70) 

a 2 nd order differential equation for £ K , an equation very much resembling New- 
ton's 2 nd law. 

Remark Note that if the electric dipole moment had the form D k E = q£, R , 
i.e., if the center of charge was the same as the center of mass, the last term 
would be exactly the classical radiation reaction force [14] . 

Returning to the Bondi linear momentum, Eq. ([6Tj) . after simplification and 
the use of Eq.l[6"4"|). becomes a dynamical expression for the momentum: 



P k = M B & - \c-\D%> + 2 -c-\KD%)'e ijk . (71) 

The dynamical or flux expression for momentum loss is obtained by taking the 
prime derivative of Eq. ([7T|) and eliminating the £ fl " by using Eq. (j70|) . i.e., 

P k ' = \c-^K^- (72) 



which is just a different form of Eq. ([70|) , 

4.4 Physical Interpretation 

We thus see that to second order in the perturbation, the imaginary part of the 
center of mass world-line leads to exactly the same equation for the radiated 
angular momentum as in the Schwarzschild case but, now, with a first-order 
(smaller) mass we do obtain recoil. The (real) center of mass (or the linear 
momentum) satisfies a 2 nd order evolution equation that is (similar to) Newton's 
2 nd law with a recoil force and a radiation reaction force. In some sense we see 
that general relativity contains Newton's 2 nd law of motion. 

In both perturbation schemes there is an anomalous contribution to the total 
angular momentum, namely in Eqs. (|49|) and (|68|l we see the term proportional 
to the rate of change of the magnetic dipole moment, i.e., D' M . This can be 
considered as a prediction though how to measure it is not clear. 

The two perturbations schemes thus lead to different physical consequences. 

5 Perturbed Metric 

For completeness, we display, to second order, the perturbed spin-coefficients 
and metric for the Schwarzschild background calculation. This metric, given 
in Bondi coordinates, is obtained by integrating the Bianchi identities, then 
integrating the spin coefficient and metric equations and finally constructing 
the metricp]. To begin, we display the leading-order radial behavior for the 
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spin coefficients; the more complicated terms are bundled together and given in 
the Appendix. 
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(73) 
(74) 
(75) 
(76) 
(77) 
(78) 
(79) 
(80) 
(81) 
(82) 



Non-radial integrating factors in these expressions can be determined from 
the relations: 



a" = -?-=£, 



(83) 
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It is useful to define an auxiliary bundled expression, 21, which enters into 
the metric expressions: 



2l = 3y»;+g^; fc g 2 k / y/sjD* + d*) | i^-^ yo 

^2 I ~.2 oJl / 1* 
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The full covariant metric to second order in the Schwarzschild pertur- 
bation, takes the form 

5oo 1 ffoc ffoc 
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5oc 5 CC ^ 
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(87) 



16 



2cr r 2a°r , s 

9C( = --j5T> 9ll = —pT> ( 88 ) 



3CC = -^T- (89) 

and 

p = (i + CO- 

The corresponding metric for the perturbation off of the Minkowski back- 
ground is easily obtained from these results. 



6 Discussion 

We begin with a mea culpa. In this work we started with two given vacuum 
metrics, the Schwarzschild and Minkowski metrics and then 'drove' or perturbed 
them both with a 'small' Maxwell field. In one case there was a 'large' mass 
while in the second case the mass was 'small'. There was no systematic attempt, 
by comparisons, to define small or large quantities; there was no small param- 
eter for a series expansion. By 'small' we meant it to be as small as needed 
for a physical effect; or as large as was needed for a different effect. The idea 
behind the two perturbation schemes was purely heuristic. We wanted to see, 
in a rough sense, within the context of the Einstein equations] 10$, what would 
be the physical responses of the recently defined (complex) center of mass of 
massive systems (with different size masses) to an electromagnetic perturbation. 
In addition to the totally expected result that the large mass object did not ex- 
perience any recoil while the small object did experience recoil, we obtained a 
variety of concomitant physical results. The imaginary part of the (complex) 
center of mass world-line could be identified with the electrodynamic induced 
internal spin angular momentum that satisfies the classical (angular momen- 
tum) conservation law. There was an electromagnetically induced gravitational 
radiation as well as linear momentum loss. We obtained explicit equations of 
motion for the center of mass that have the form of Newton's 2 nd law. 

The results presented here, in the two different models, are in total agree- 
ment with well-known classical Newtonian and electrodynamics effects. They 
however go beyond these known results. They give a confirmation of the physi- 
cal ideas developed purely in the context gravitational theory (general relativity) 
for the identification and the associated dynamical response of certain geomet- 
ric structures that arise naturally. The most important of them is the special 
class of null geodesic congruences, the shear-free (or asymptotically shear-free) 
null geodesic congruences. It is their very existence that yields or provides the 
complex Minkowski world-lines that become our complex center-of mass which 
then yield our angular momentum expressions. 

As a final comment we point out that the perturbed space-times described 
here contain, as special cases, both the Kerr and the charged Kerr metrics as 
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perturbations off Schwarzschild. They occur when the electromagnetic dipoles 
are 'shut off', i.e., when D l = 0. 
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9 Appendix 

9.1 Radial Behavior of the Weyl Tensor 

Below are the full expressions for the bundle terms (Ai) given in Eq. [Mil 

A 1 = --i2 q D*Y 1 \ + 2V2D*&'( — e ljk Y 1 \ + -Y 2 \ ] )} (90) 
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^3 = -^r + ^r-^-Yu 1 (92) 

„, , (1 + 7V2)D 1 ' + 6V2D l/ (IID' + UD 1 ) 2V2GM S D\„_ 1 
+2kq{ 3^ + I&3 ^—) Y u 

(34 - V2){D l &)' - D { D j ' hV2D l D3 56GM S DW^ _ x 

+ ( 9^ + ~~ 4^ 5^ )tijkYlk 



^jyJ'jyi 20D l 'Di' + y/2{D i D i ')' , 
"374 



,y au u j z,uu u j -t- v ^\i-> i-> ) \tt-i 



•il^/2D 1 '& + (3 - 49\/2)(.D l i>)' 23D l D^> x 
( 18^ + 12r6 )y2i J 

2V2Gfc V2M s D l D3' 16M S D 1 D^ x 

^ ( ~& + 5^7 ) 2ij 



^ - -— + ^r--&r _ ( 93 ) 

3V2D l "Di' + 2{D t "&) 1 (16\/2 - 72)(L> i 'i>)' 
~ I ^ ^4 

128D l 'Di' (24-107V2)(D l Diy -25V2D l 'Di A1D 1 & 

9T 3 + 36^ + 15r 6 

2A0GM s D ll & 38V2GM s D l & 2 

+ l 5c 2 r 6 ^7 ) «j ' 

9.2 The Spin Coefficients 

The following are the full expressions for the &i given in Section V for the spin 
coefficients in the Schwarzschild background perturbation: 
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& P = -^ Y ^ + lk (^^-^^^ kY ^ ( 96 ) 
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